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The pseudogap phase is considered as a new state of matter in the phase string model of the 
doped Mott insulator, which is composed of two distinct regimes known as the upper and lower 
pseudogap phases, respectively. The former corresponds to the formation of spin singlet pairing 
whose magnetic characterizations have been recently studied [Phys. Rev. B 72, 104520 (2005)]. The 
latter as a low-temperature regime of the pseudogap phase is systematically explored in this work, 
which is characterized by the formation of the Cooper pair amplitude and described by a generalized 
Gingzburg-Landau theory. Elementary excitation in this phase is a charge-neutral object carrying 
spin-1/2 and locking with a supercurrent vortex, known as spinon-vortex composite. Such a lower 
pseudogap phase can be regarded as a vortex liquid state due to the presence of free spinon-vortices. 
Here thermally excited spinon-vortices destroy the phase coherence and are responsible for nontrivial 
Nernst effect and diamagnetism. The transport entropy and core energy associated with a spinon- 
vortex are determined by the spin degrees of freedom. Such a spontaneous vortex liquid phase can 
be also considered as a spin liquid with a finite correlation length and gapped S = 1/2 excitations, 
where a resonancelike non-propagating spin mode emerges at the antiferromagnetic wavevector (vr, tt) 
with a doping-dependent characteristic energy. The superconducting phase is closely related to the 
lower pseudogap phase by a topological transition with spinon-vortices and -antivortices forming 
bound pairs and the emergence of fermionic quasiparticles as holon-spinon-vortex bound objects. 
A quantitative phase diagram in the parameter space of doping, temperature, and magnetic field is 
determined. Comparisons with experiments are also made. 

I. INTRODUCTION 

Pseudogap phenomenon^'^ in the underdoped cuprates is widely considered to be closely related to the high- Tc prob- 
lem. Any sensible microscopic theory for the high- Tc superconductivity is expected to include the pseudogap as an 
integral part of the theory concerning the high-temperature/energy and short-distance physics. Due to the unconven- 
tional properties manifested<li2, in transport, spin dynamics, optical and single-particle spectroscopies, thermodynamic 
properties, etc., experimental measurements in the pseudogap regime have imposed much stronger constraints on a 
potential high- Tc theory than just those in the superconducting state. 

Experimentally two distinctive pseudogap regions have been observedi*^*^*^. The high-temperature (T) regime, 
called the upper pseudogap phase (UPP), is marked by the suppression of uniform magnetic susceptibilitji2i& and the 
deviation of dc resistivity from the linear-T behavioB^iSiiS, both of which happen below a characteristic temperature 
Tq. The temperature Tq is comparable to the superexchange coupling strength (~ 1,000 K) near the half- filling and 
decreases monotonically with dopin g'^i^i^d" . Near the optimal doping, Tq is reduced to the same order of magnitude 
as the superconducting transition temperature Tc. The single-particle pseudogap structure around the wavevector 
(±77, 0) and (0, ±7r) observed in the ARPES experimentsii also shows a similar doping dependence. A low-T pseu- 
dogap regime, called the lower pseudogap phase (LPP), corresponds to a crossover regime between the UPP and 
superconducting (SC) phase, where the low-energy antiferromagnetic (AF) fluctuations get suppressed as shown in 
^^Cu NMR experiments^ii^, accompanied by the emergence of a "high-energy" ("41 meV") resonancelike peak in the 
dynamic spin susceptibility function around the AF wavevector (tt, tt) as revealed by the neutron measurements^^. 

One of the most significant experimental observations in the LPP region is the discoverjiMiiSii24iLiS4i2i2£ of the 
large Nernst signal and residual diamagnetism which strongly support the presence of vortex excitations2i in a wide 
temperature range which can be several times higher than Tc at low doping. Generally it is believed that the 
superconducting phase transition in underdoped materials is due to the phase disordering effect^^ by vortices as one 
approaches Tc from below. There have been a number of theoretical proposals^SiSiSiiSSiSSiSi to explain the existence 
of topological vortices above Tc in the underdoped regime. However, the nature of these vortices above Tc can vary 
distinctly from the more conventional Kosterlitz-Thouless (KT) vortice3S2i2SiS& to exotic vortices with a spin-1/2 
(spinon) trapped at each coreS^. 

It would be rather difficult for one to judge the relevance of a particular microscopic theory simply in terms of some 
detailed account for individual experiments, considering huge amounts of existing data available for interpretation and 
the lack of consensus on the underlying mechanism after almost two decades' extensive and intensive investigations. 
It is thus particularly meaningful for a potential microscopic theory of high- Tc to be worked out extensively and self- 
consistently to provide a comprehensive "road map" . One may then use the obtained specific "map" and organizing 
principles to accommodate and analyze the experimental results with priority and hierarchy, and to further predict 
new experiments. Only in this way may one effectively justify or falsify a particular theoretical approach based on a 
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FIG. 1: The global phase diagram of the phase string model at low doping—. The high-temperature phase is called the 
upper pseudogap phase (UPP) which is characterized by the bosonic RVB order parameter A" and the temperature Tq. The 
antiferromagnetic (AF) long-range order state, the lower pseudogap phase (LPP) which is also called— the spontaneous vortex 
phase (SVP) with a characteristic temperature r„, and the d-wave superconducting phase (SC) below Tc, all occur as some forms 
of low-T instabilities in the UPP. For example, the AF state corresponds to the spinon condensation on the RVB background, 
while the LPP and SC states involve the holon condensation. The detailed phase diagram is determined by the unique mutual 
duality gauge structure—*— which quantitatively governs the intrinsic competition between the charge and spin degrees of 
freedom. The shaded regime of the LPP will be the main focus of this work. 



broad set of experimental measurements. 

The "road map" of the phase string model^'* is illustrated by a global phase diagram at low doping in Fig. 1. 
This model is an effective theory proposed for the cuprate superconductors based on the t ~ J model, with a careful 
treatment of the spin correlation at various ranges and with self- consistently incorporating the mutual interplay 
between the spin and charge degrees of fj-eedomS^iSSi^fliSiiSi^^. Such a model describes a doped Mott insulator with 
spins forming singlet pairing, known as the bosonic resonating-valence-bond (RVB) pairing, in the UPP. The magnetic 
characterizations of this phase have been carefully examined in Ref.^^ recently. At half-filling, without the presence of 
holes, the short-range AF correlations in the UPP will continuously grow with reducing temperature and eventually 
becomes AF long-range ordered (AFLRO) in the ground state (a finite Neel temperature will be the consequence of 
the interlayer coupling). On the other hand, as doped holes moving in the spin background always generates the 
so-called phase string effect^^, the spin dynamics"^'"' will be drastically reshaped beyond some doping concentrationSi, 
in order for holes to further gain kinetic energy^^ at low-T. In such a finite doping regime, one enters the LPP [also 
known as the spontaneous vortex phase (SVP)^] at a low temperature T„ and finally ends up in the SC state2L2& 
below Tc without encountering an AFLRO. 

The above global phase diagram with the UPP as the high-T "unstable fixed point" state of the doped Mott 
insulator is spiritually very similar to the proposals made by Anderson^LSS recently, where the LPP and SC phases 
are considered as low-energy /temperature instabilities with the spin-charge "locking" evolving progressively. Although 
two approaches look different mathematically, they still share some important features like the RVB spin pairing in 
the UPP, vortex liquid state in the LPP, and vortex-antivortex binding in the SC, as well as the kinetic-energy-driven 
and progressive charge-spin recombination in the latter regimes. The main distinction between the two approaches 
perhaps lies in how the UPP is described as the high-energy/temperature starting point. In a framework based on 
the Gutzwiller-projected fermionic RVB (f-RVB) mean-field stat o^i"^° , the short-range AF correlations are normally 
intrinsically weak, before the Gutzwiller projection, as compared to those in the bosonic RVB state, such that the 
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kinetic energy of doped holes at the mean-field level is not as severally frustrated. The advantage of this approach 
is that the coherent fermion nodal quasiparticle appears directly in the formalism. By contrast, in the bosonic RVB 
state, the short-range AF correlations are intrinsically associated with the bosonic RVB pairing, which remain strong 
even before the Gutzwiller projection is made or in other words, at the mean-field level, such that the motion of 
doped holes always gets severely frustrated in the UPP. Consequently, with decreasing temperature, the tendency 
for gaining the kinetic energy becomes the critical driving force to realize the LPP and SC at a finite doping. Here 
the fermion nodal quasiparticle emerges as a coherent one only in the SC state through a more delicate spin-charge 
binding process*^. It is also noted that the gauge structure describing the fluctuations around the mean-field state, 
which is either 11(1)^24^, SU(2)^^ or even Z^^, arising from the no double occupancy constraint in the slave-boson 
approaches24£ based on the f-RVB description, is rather different from the present mutual duality (mutual Chern- 
SimonsS) gauge structure arising from the phase string effect in the t — J modelS. The latter is crucially responsible 
for the global phase diagram shown in Fig. 1, occurring at proper temperatures and doping regions as compared to 
the experiments. 

In this paper, we will present a systematic and detailed description of the LPP (the shaded region in Fig. 1) based 
on the phase string model in the context of the global phase diagram involving the UPP, SC, and AF phases. The 
basic picture is that the pseudogap phase is a new state of matter, which is a doped Mott insulator with spin-charge 
separation. In particular, the LPP is characterized by the holon condensation, described by a generalized Gingzburg- 
Landau (GL) theorySSiSi. The elementary excitation identifled in this phase is a bosonic charge-neutral object carrying 
spin- 1/2 and vortexlike currents known as a spinon- vortex. The LPP is composed of thermally excited spinon- vortices 
which destroy the SC phase coherence, while the Cooper pair amplitude still remains. Interesting consequences of this 
phase will be explored, including the Nernst effect and diamagnetism. How the LPP is closely related to the SC phase 
by a topological transition involving vortex-antivortex binding or "spinon confinement" , including the emergence of 
coherent fermionic quasiparticles in the SC phased, will be also discussed. 

The remainder of the paper is organized as follows. In Sec. II, the theoretical framework of the phase string model 
is presented and the elementary equations describing the LPP are given. In Sec. HI, the physical properties of the 
LPP are explored, including the spin liquid behavior, elementary excitation, phase diagram, the Nernst effect and 
diamagnetism, as well as the topological transition to the SC state. Finally, a summary and conclusions are presented 
in Sec. IV. 



II. ELEMENTARY EQUATIONS 
A. Phase string model 

The phase string model is an effective low-energy theory, obtained^2i2^ based on the phase string formulation of 
the two-dimensional (2D) t — J model, in which the electron annihilation operator is decomposed in the following 
bosonization form^ 

= /il&,.e'^- (1) 

where both holon h}^ and spinon bi^ operators are bosonic fields, satisfying the no-double-occupancy constraint 
n-i + — 1 with — h\hi as the holon number and — b\^biu as the spinon number. The phase factor 

e*^'", which ensures the fermionic statistics of cia, is defined by e*^'" = e'^^^" " (ct)^'' (— ct)*, where Nh is the total 
holon number operator, a = ±1, and (—1)* = ±1 (i € even/odd site) is a staggered sign factor. The phase Q^^"^^ 
= \ [$1 — <i>° — 0-$^] is a nonlocal operator with 



(2) 



*° = T.l(^^) Im In [z, - zi), and 



= ^ Im In {z, - zi)4 , (3) 

where zi — xi + iyi is the complex coordinate for a lattice site I. 

The effective phase string theory is a bosonic RVB theory, in which the spin degrees of freedom are described by 
the following spinon Hamiltonian: 
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Hs ^ -Js E [(^'^^-)c^t.+H.c; 

<ij>a 

-AE^U'^-^) ■ (4) 

Here the effective coupling constant Js — J{1 ~ g)A'' /2, where J is the superexchange constant in the original t ~ J 
model, g ~ 45 with d denoting the doping concentration stands for a renormalization effect to the superexchange 
energy and A'* is the bosonic RVB order parameter defined by 

with the sites i and j belonging to two nearest neighbors (NN). The Lagrangian multipher A enforces the constraint 
Ea{blb^a) = l-S. 

The charge degree of freedom is governed by the holon Hamiltonian 

Hk = -thYl i/^'^^'"^'^) h\hj + H.c. (6) 
<»i> 

where th ^ t {t is the hopping integral in the original t — J model) is the renormalized hopping integral for holons, 
which directly see the external electromagnetic vector potential Al^ with charge +e. 

A peculiar feature of the phase string model is that the spinous and holons see two different gauge fields, A^^ and 
Alj , with their strengths constrained to the densities of holons and spinons, respectively. Namely, A'^'j and Afj satisfy 
the following gauge-independent topological constraint 

/dr- A'' =7riV''(c) (7) 



j( dr • A'* = TT [iVj- (c) - A^f (c)] (8) 

where A'' and A* are defined by A^j = Yij -A^ and Afj = Yij •A'', respectively, with (|ry | = a is equal to the lattice 
constant) a spatial vector connecting two NN sites, i and j. On the right-hand side (rhs) of the above expressions, 
N^{c) and i^{c) denote the numbers of holons and spinons of spin "f or |, respectively, which are enclosed by an 
arbitrary counter-clockwise-oriented path c. Such a topological gauge structure entangles the otherwise decoupled 
spin and charge degrees of freedom and is called mutual duality since a spinon sees a holon as a tt flux tube and vice 
versa, which has been also formulated in the path-integral formalism as a mutual Chern-Simons gauge theory— with 
both time-reversal and spin rotational symmetries retained. 

The phase string model with the mutual duality gauge structure can lead to a global phase diagram at low-doping as 
summarized^^in Fig. ^ It been shown'^'* that as the temperature is reduced below the characteristic temperature 
To of the UPP, the spins in the background start to form singlet pairs (bosonic RVB pairs) and condense, which is 
characterized by the RVB order parameter A* . The short-range AF correlations simultaneously begin to grow below 
To, in contrast to much weaker spin correlations above Tq. 

Several low-temperature phases, including the so-called LPP, which is also known as the SVP, the d-wave super- 
conducting phase, and the AF phase at low doping, are all nested below the UPP as shown in Fig. 1 with essentially 
the same symmetry in the spin sector as characterized by A*. The competitions between these low-temperature 
phases are basically decided by the above mutual-duality (mutual Chern-Simons) gauge structure, rather than more 
conventional competitions between different order parameters. 

In particular, the LPP corresponds to the occurrence of holon Bose-condensation i^hl^ 7^ without superconducting 

phase coherence, whose region is shaded in Fig. ^and will be the main focus of the following study in this paper. 



B. Superconducting order parameter 

The LPP is the holon condensed phase on a background where spins form the singlet pairing. In order to understand 
the corresponding physical characteristics, it is instructive to examine the superconducting order parameter A^^^ = 

O- {CiaCj-cr) bclow. 



(a) 



(b) 



FIG. 2: (a) An isolated spinon is always associated with a 27r vortex in the phase of the superconducting order parameter in 
the phase string model, (b) A pair of vortex-antivortex with spinous located at the cores. 

Based on the phase string formulation [Eq. (^], the spinon RVB paring [Eq. ISJ], and the holon condensation 
condition (hl^ Oj the superconducting order parameter defined at the NN sites can be expressed aaS-'^ ^sc = 

description we shall treat the charge condensate in terms of a slowly varying continuous field, (hf^ ~^ ''l^hi^i)^ is 

convenient to rewrite ^'sc in the continuum limit (without considering the d-wave symmetry of the relative coordinate 
for simplicity) as follows 

vl/scoc A0(e'^*''M) (9) 

where the amplitude 

A° = A^(VD' (10) 

and the phase $^(ri) = $f is given by 

$''(r) = J d?Y' Im In [z ~ z'\ {n\{v') ~ n\{v')\ (11) 

withn^(rO = nt/a^ 

Therefore, in the phase string model the superconducting order parameter ^'sc has an intrinsic composite structure 
of amplitude and phase components in Eq. (Q. Here the Cooper pair amplitude A*' is composed of the RVB pairing 
A* and holon condensate which is finite at T <Ty [note that generally Ty < Tq, a,s tph 7^ is always underpinned 
by the spin singlet paring in the present work (c/. Sec. Ill C)]. Thus the LPP can be physically regarded as the 
formation of the Cooper pair amplitude A'^ 0. 

However, the onset of A" does not lead to superconductivity immediately. From Eq. (|ll|l . it is clear that 
describes 2tt phase vortices whose cores are centered at spinous: i.e., — s- (f'* ± 27r or ^sc ^'sce^^'^* each time as 
the coordinate r continuously winds around a spinon once according to Eq. (|ll|l . In other words, a spinon is always 
associated with a 27r vortex in known as a spinon-vortex composite which is schematically illustrated in Fig. 

|2Ia). A spinon-vortex and -antivortex pair of a finite range of separation will result in the cancellation of the phase 
at a large length scale as shown in Fig. Ofb). Consequently, if there exist excited spinous freed from the RVB 
pairing, then is generally disordered such that (^e**"^''^) = and "ifsc — 0, which defines the LPP at Tc <T < Ty. 
Here the presence of free 2tt vortices in the LPP is a direct manifestation of the electron "fractionalization" , i.e., the 
existence of neutral S = 1/2 spinous in the underlying phase. In the opposite, the SC phase is determined by the 
phase coherence (e** *•'"■') ^0 which corresponds to the spinon pair confinement with no single spinon excitations 
allowed in the bulk of the sample in the Meissner phase. Thus, the LPP and the superconducting phase are very 
closely related, distinct mainly by a topological transition in the long-distance limit which is to be further discussed 
in Sec. HI later. 
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Finally, we point out that a spin 1/2 bound to the core of a 27r vortex of the superconducting order parameter is a 
very general and natural consequence of doped Mott physics. Namely, with the depleting of superfluid (hole) density at 
the core center, due to the no double occupancy constraint, a spinon trapped inside is energetically most competitive. 
Generally there could be also two other ways to create topological excitations in the superconducting order parameter 
according to Eqs. (jSJ and H10|l . which correspond to a 27r- vortex singularity in the holon condensation order parameter 
■0/1 (r) or the RVB order parameter A*, respectively. However, the former vortex defined as "0/1 (r) ■)/)/i(r)e*^"^'"^^^^'') 
leads to a An vortex in ^'sc, which is energetically not favored. (Nevertheless, such Att vortices do play a role when 
they are bound to the spinon vortices, as will be shown in Sec. II C.) The latter one does provide a 27r vorticity in 
^SC, but it costs an expensive core energy due to the suppression of the RVB pair order parameter. A'' 0, at the 
vortex core. By comparison. A** can still remain finite in the core of a spinon-vortex shown above, where the average 
^gj<i>,(r)^ is suppressed to zero within a scale defined by an intrinsic length scale of the unpaired spinon trapped inside. 
As to be shown later, the existence of spinon-vortices in the present theory provides a natural origin for "cheap" 
vortices in the LPP, with the core energy related to a characteristic spin gap scale. 



C. Gingzburg-Landau equation 

According to the above definition, the LPP represents a state of matter characterized by an order parameter field 
'0ft,(r), on top of the bosonic RVB condensate, whose microscopic definition is given by 

Mr) = {h{r)) (12) 

where h{r) is the continuum version of the bosonic holon matter field experiencing the Bose condensation. Different 
from an ordinary thermodynamic phase like the BCS superconducting state, in which the order parameter is expressed 
in terms of the electron operators {i.e., a pair of the electron creation or annihilation operators), the holon operator 
h{r) cannot be simply expressed based on a local combination of the electron operators. The corresponding off 
diagonal long range order (ODLRO) as sugg ested by Eq. (O: {h^r)h{r')) ■0*(r)0,,(r') at |r-r'| oo, must 
be then a "hidden" ODLRO, meaning that it involves a nonlocal, infinite-body order parameter if to be expressed in 
terms of the electron operators. More profoundly, since iphi^) will couple to an internal gauge field whose symmetry 
cannot be spontaneously broken, the onset of Eq. (|12|l . or more precisely, its amplitude (see below), should actually 
represent a crossover rather than a true phase transition. 

In the following we introduce an effective GL description of the LPP based on the order parameter iph. Let us start 
with the continuum version of Eq. J^Jl which is given by 

Hh = — / d^r h\r) (-iV - A" - eA'f h{r) (13) 
2m/i J 

where ruh = {2tha'^)~^ and A'^ is the electromagnetic vector potential which is seen solely by the holons which carry 
charge +e. By noting that the holons here are hard-core bosons with repulsive short-range interaction, one may 
generally write down the corresponding GL free energy Ft = J (Pr fh where^ 

h = + a + ^|0,|4 + -^0* (-zV - A^ - eA-^)" 0, (14) 
2 Irah 

with denoting the "normal state" free energy density. In the absence of fields and gradients in Eq. (|14|l . the 
minimum of fh = fjl + a iV-'hl^ + f gives rise to 

at a < 0. Here p° represents a bare "superfluid" density, which is assumed to remain finite throughout LPP and is 
reduced to a renormalized ph = \tph\^ in the presence of A*. 

By minimizing F^, a nonlinear Schrodinger equation or differential GL equation satisfied by i^hi^) can be obtained 

ai^h + vM^h + (-iV - A^ - eA'f 0;, = (15) 

with the standard "supercurrent" density given by 

7 A -I- p A ^ 

J(r) = -TT- mr)^Mr) - V0^(r)V,(r)] ± 0^(r)0.(r) . (16) 

2mh mh 
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These equations are similar to an ordinary GL theory describing a charge +e Bose condensate coupled to an external 
electromagnetic field A^, except that iph is further coupled to the spin degrees of freedom through the vector potential 
A*^, whose gauge-independent definition is given in Eq. 0. It means that each isolated spin (spinon) will register as 
a ±7r flux tube in the nonlinear Schrodinger equation l|15|) to exert frustration effect on the charge condensate, which 
reflects the key influence of the spin degrees of freedom on the charge condensate in the phase string model. 

The GL equations, Eqs. H15|l and H16|l . are gauge invariant under the transformation: iph ~^ iphe^"^ and A** 
A^+Vip where ip is an arbitrary single- valued function. Especially, one may also construct a large gauge transformation 
'P^J^zVi, with 

dr ■ Vlp,. = ±27r (17) 

where a loop q encircles a ±27r vortex center i which coincides with a spinon position. This procedure is allowed 
because in the underlying microscopic theory a holon cannot stay at the sites of spinons due to the no double occupancy 
constraint and iph should thus vanish at a spinon position where a ±7r flux tube is bound to according to Eq. (jHJ (A" 
diverges there). Through such kind of large gauge transformation, the rhs. of Eq. © (if one spinon is enclosed by c) 
may be properly transformed by a minus sign as 

dr • A" = ±7r ^ j( dr • A" = TTT (18) 

which means the sign of vorticity around a spinon is no longer directly tied to the spin index a =|, [. In this way 
the theory is expressed in an explicitly spin-rotational invariant form, where the sign of the tt fluxoid carried by each 
spin excitation (spinon) can be determined by minimizing the free energy instead. 

The above large gauge transformation can be also visualized based on the superconducting order parameter 
defined in Eq. 0. Evidently, <I>^ is related to A* (V^'* — 2A'*), satisfying the gauge invariancc under the gauge 
transformation: ■0fi iphe'^''^ , A^ A^ -I- Vip, and -|- 2(p. If one makes large gauge transformations like Eq. 

((T7|l . then the vorticity signs of the 2tt vortices in can be changed such that ^^(r) and n^^{r) generally should be 
replaced by 

^'{r) ^ J d^r' Im\n[z- z'] [n'L(r') - «+(r')] , (19) 

where ri±(r) denotes excited spinon number with vorticity ±27r no longer be directly tied to the spin index (the 
definition of ± sign for the vorticity is such that -I- vortices are commensurate with the magnetic field applied along 
the perpendicular direction of the 2D plane, cf. Sec. Ill B). In order words, although (V'^)^ only contributes Att 
vortices to ^sc, it can also change the sign of 2tt vortices in if the vortex cores of these Att vortices properly 
coincide with the spinons (note that the higher vorticities like Qtt per spinon as the result of the above transformation 
will cost too much energy and are thus not considered here). 

In summary, the LPP phase is described by a generalized GL description of the holon condensate coupled with 
unbound spinon- vortices, the dynamics of which is governed by the spinon Hamiltonian In the next section, the 
physical properties of the LPP phase will be explored based on such an effective description. 

III. PHYSICAL PROPERTIES OF LOWER PSEUDOGAP PHASE 

Let us first focus on the spin degrees of freedom below, which will play an important role in the LPP. 

A. Spin degrees of freedom: A spin liquid 

Let us recall that the spinon mean-field Hamiltonian |0J can be diagonalized by the following Bogoliubov 
transformationSii 

ha = ^ W^a (r*) Umlma " VnilL-a (20) 
m 

with u,n = l/\/2(A/£'„i + 1)^/^ and Vm — i/V^{^/E„i — l)^/^sgn(^„i). The diagonalized Hamiltonian is written as 

Hs = Y^ E^iclLalnia + const. (21) 
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where the spinon spectrum Em = \/X^ — S,^^ ^^d is the eigenvalue of a tight-binding equation 

^mWrncr{ri) ^ ~Js ^ e~'-'^^'i"Wma{r j) (22) 

Here the gauge field A'^-j is constrained to the hole number according to Eq. , which incorporates the most nontrivial 
doping effect on the spin degrees of freedom in the phase string model. In the LPP when holons are Bose condensed, 
A'^j may be treated as describing a uniform flux of strength Stt per plaquette. 

A spinon excitation created by 7,J^g. with an eigen-energy E„i has no direct gauge coupling with the external 
electromagnetic field A'^ in the mean-field spinon Hamiltonian Hs [Eq. H21(l ] and is thus a charge- neutral, 5' = 1/2 
object. Here a magnetic field can still couple to spinous via the Zeeman energy, — /xb-B X^icr ^"-io-- Adding this term 
to Hs will lead to a modified spinon spectrum in Eq. 121|l 



Ema — Em — CMB-B (23) 

In obtaining Eq. H23|) one need to note that 



ma 



with using — = 1 and ^ |wmcr(ri)|^ = 1, where the excited spinon number rt^^ is defined by 

i 

'"^Ina = llncrlma- (25) 

In the LPP, by treating A^^ as a vector potential for a uniform flux perpendicular to the 2D plane with the field 
strength = the spinon wave function Wma (r) determined by Eq. (|22|l is characterized by a wavepacket of a 
typical size of the cyclotron length scale Oc — l/^B^, i.e., 

ac = (26) 

V TTO 

The RVB pairing of spins described by Hg in Eq. (|21|l has the following mean-fleld form in the ground state^^ 

\mB)=Ce^v{Y.W,,h\^h]}j |0) 

where the RVB amplitude \Wij \ — ^w^a (r^^mcr (i"j) behaves as follows 

at |ri — Fj l 3> a, for ij belonging to different sublattices. Here the characteristic size for the RVB pair is given hy^ 

Cs = a/5 (27) 

which also decides the equal-time spin-spin correlation length. Note that diverges at small 5 where an AFLRO 
emerges. Namely, the spin sector is described by a spin liquid at finite doping, which can be continuously connected 
to the AFLRO state as the doping concentration is reduced to zero. 

Thus, although there exist a lot of spins at low doping {\ — 6 per site), the majority of them form short-range RVB 
pairs of a typical size in the holon condensation phase. According to Eq. then the contribution of these RVB 
paired spins to A* will be effectively cancelled out, and the charge condensate described by iph is essentially decoupled 
from the RVB background at a scale larger than ^s. This is also self-consistent with the fact that is compatible 
with the average holon-holon distance which sets a minimal length scale for the holon condensate. Only the excited 
spinous will then effectively contribute to A* in the GL equation by noting that 

E^("-) = E^l^-(^')l' (28) 
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FIG. 3: Spin liquid behavior in the LPP ai 5 — 0.125 as compared to the half-filling case with an AFLRO. (a) Dynamic spin 
susceptibility at momentum Qaf = (i",""). Inset: the evolution of the resonancelike peak energy Eg as a function of doping 
concentration, (b) Uniform spin susceptibility versus the temperature. 



based on the Bogoliubov transformation H20|l . In obtaining Eq. H28(l . one has used the fact that |wmcr(r)| is 
independent of a. Here the length scale determined by \'Wma{''^i)\ , «-e., the "cyclotron" length a^, in Eq. 
sets a natural spinon size which is within the "coarse graining" length as Oc = ^s/\/2 ^ ^s- Thus an excited spinon 
can still maintain its quantum integrity in the GL equation description. 

Eq. H21|l describes the spinon excitations as neutral spin-1/2 excitations in the RVB background, labelled by m 
and spin index a. Fig. OJa) shows the corresponding dynamic spin susceptibility function at the AF wavevector 
Qaf = (tt, tt) for 6 = and 6 = 0.125, respectively, in which a resonancelike peak naturally emerges at finite doping. 
In the inset, such a resonance energy Eg as a function of doping is presented (the bare superexchange coupling J 
used in Ref.-^ should be replaced by a renormalized one Jeff in the expression of Jg as discussed in Ref.^"*). Note that 
more precisely, the competing AF state can drive Eg vanish at Xc ~ 0.04 in a fashion of Eg ~ ^/S — XcJ,^^ instead of 
Eg ^ 6J a.t 6 —> a.s shown in Fig. Ofa), which is beyond the mean- field description presented here. In Fig. Ofb), 
the uniform spin susceptibility'^^ is also shown for both 6 = and 6 = 0.125, and the later one clearly exhibits a 
pseudogap behavior at low temperature. These typical magnetic properties evolving from the undoped AF phase to 
the LPP are quite consistent with the experimental measurements in the high- Tc cuprates. 



B. Spin-charge entanglement in the LPP: Spinon- vortex composite 

As discussed above, the spin sector in the LPP can be characterized as a spin liquid state with charge neutral, 
spin 1/2 spinous as its elementary excitations. Here one finds essentially the same symmetry for the spin degrees 
of freedom in both the UPP and LPP, which are described by the bosonic RVB order parameter A*. What really 
differentiates the UPP and LPP lies in the charge sector by the additional holon condensation -0/1 ^ in the LPP. 
As it turns out, this distinction in the charge degree of freedom then feedbacks to the spin sector to affect the spin 
excitations, resulting in a spin gap behavior in the LPP as represented in Fig. |3in the previous section. 

In the following we shall discuss a more striking feature due to the mutual entanglement between the spin and 
charge degrees of freedom based on the GL equations. Namely, a spinon excitation will always induce a current vortex 
in the LPP to form a spinon-current-vortex composite, which is consistent with the 27r vortices previously identified 
in the phase of the superconducting order parameter in Sec. II B [Fig. 2(a)]. 

By writing 



(29) 
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FIG. 4; A neutral spin-1/2 bosonic object, spinon, is always locking with a current vortex as an elementary excitation in the 
LPP, known as a spinon- vortex composite. The vortex core size is decided by the quantum length scale of the spinon, Oc, with 
its vorticity sign independent of the spin index. 



a London equation for the supercurrent based on Eq. H16|l (neglecting the spatial dependence of pu) is given by 

J(r) = [Vc^h - A' - eA^] (30) 
rrih 

Since each unpaired spinon will contribute to dr ■ A* = ±7r in terms Eq. (jH)) if the loop c encloses such a spinon of 
a typical radius of Oc, a minimal supercurrent vortex centered around it is then given by 



/dr • J(r) = ±7r-^ (31) 
Jc ruh 



according to Eq. H30I) . This spinon- vortex composite is illustrated in Fig. 01 where the "cyclotron" length Uc serves 
as a natural radius for the core size of a spinon- vortex. 

As discussed in Sec. II C, a large gauge transformation in (j)^ can result in a binding of a 2tt vortex to a spinon to 
change the sign of vorticity it carries. Consequently, a spinon excitation is always a composite of a neutral spin-1/2 
locking with a charge current vortex in the LPP, but the vorticity and spin index are not necessarily the same. Such 
a spinon-vortex binding is a key consequence of holon condensation and play an essential role in understanding the 
LPP physics. Due to the presence of unconfined spinon-vortices, the LPP is a phase-frustrated BEC of holon and is 
more precisely characterized by the finite amplitude of iph , ^/ph- 

By requiring that 



j(dr • J(r)^ = 



and omitting independent 27r vortices in V^/i (which are presumably important only at temperatures close to T^,), 
the following "neutrality" condition is then obtained 

(fdr-A^ = - I dr ■ {\7(j)h- A") 
Jc e 

= MN+{c)-N^{c)] (32) 

where = hcjle and N±{c) denotes the number of spinon-vortices ( ant i vortices) enclosed in an arbitrary loop c. It 
shows that a perpendicular magnetic field will generally polarize the numbers of ±7r spinon-vortices (the sign ± is 
defined according to Eq. (|32|l 'l in the LPP. 

Thus, spinon excitations in the LPP are no longer "free" ones solely determined by Hg in Eq. (|21|) . They will 
actually acquire a long-range logarithmic interaction due to the spinon-vortex effect. Indeed, substituting Eq. H29|) 
into Eq. H14|) and carrying out the area integration, one finally obtains 

Fh = J J d^ridh^ [l]^a«a(ri)] ^^(ria) [Y. f^'^li''^)] + const. (33) 
in which a, /? = ± refer to the signs of vorticities carried by spinous and 

Vir,,)^-^ln\^^ (34) 



with Tc a. Nevertheless, in the LPP where a substantial number of spinon-vortices is present, the long-range force 
in Eq. (|33|l will be screened such that spinon-vortices are not bound together, as opposed to the lower-temperature 



11 




FIG. 5: The characteristic temperature and magnetic field scales which determine the phase diagram of the upper and lower 
pseudogap and superconducting phases based on the phase string model. Here To and HpQ = Hpg{T — 0) for the UPP, Tv 
and H° = H^{T = 0) for the LPP, and Tc for the SC phase. 



superconducting regime where all excited spinous form vortices-antivortices bound pairs (c/. Sec. Ill F). Namely the 
LPP remains an electron fractionalized state before it undergoes the superconducting phase transition at Tc- 

However, in the LPP, the short-range interaction in Eq. remains unscreened. Each excited spinon at a state 
labelled by m has a spatial extension determined by \wm (r) | as mentioned before. With involving the charge current, a 
symmetric profile of |wm(r)| in the symmetric gauge is energetically most favorable. Then a pair of spinon excitations 
at the same m with opposite vorticities have essentially no contribution from Eq. H33|) due to the cancellation, but 
one with the same vorticities will cost a big repulsive core energy ~ in Fh- Such a short-range interaction effect 
will be further considered later when we discuss the diamagnetism in "tfee LPP. 



C. Phase diagram 

The phase diagram of the doped Mott insulator can be drawn by combining the description of the LPP phase here 
and the previous work on the UPP phase"^'*. 

The upper pseudogap phase. — The phase string model holds for the whole bosonic RVB regime underpinned by 
A* 7^ 0, which defines the UPP as illustrated in Fig. 1. The characteristic temperature Tq and magnetic field Hp^ 
(at T — 0) for the UPP at the boundary A'* = have been previously determined in Ref.'^'*, which are consistent 
with the experimental measurementaiSiSiiS. Such temperature and magnetic field scales are replotted as functions of 
doping with xrvb = 0.25 in Fig. O (filled triangles). 

The lower pseudogap phase. — Now let us consider the phase boundary of the LPP, which is set by vanishing -0^ at 
Ty. Note that without A** this would be a simple 2D hard-core boson problem according to Eq. Ijl^f) . and the KT 
temperature for the holon condensation is given by 

Tkt = nS {2a^mh)~^ (35) 

which has been previously used^^ as an upper bound temperature for the LPP at low doping, shown in Fig. 5 by a 
dotted line (with th = 3J). 

However, the frustration effect of A* on the holon condensation can be very important. The spinon- vortex density is 
determined by the density (per site) of spinon Bogoliubov quasiparticles as n^, = J^ma ('^mo-) Due to the opening 
up of a spin gap Eg in the holon condensation phase {i.e., the LPP), n„ is exponentially small for T <C Eg. With 
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FIG. 6; Two holons (open circles) exchanging positions can pick up a minus sign if an unpaired spinon (grey circle) is enclosed 
in between, which carries a ±7r flux tube due to the mutual duality gauge structure in the phase string model. So the phase 
coherence of bosonic holons can be effectively destroyed to prevent the Bose condensation if there is a sufficient number of free 
spinon excitations in the singlet spin background. 

increasing temperature, n„ will monotonically increase until reaching the maximal number rt^^'^ = 1 — 5 at T = To 
where all the RVB pairs break up. 

Since each free spinon carries a tt fluxoid as perceived by the holons, the quantum phase coherence among bosonic 
holons will be violently interrupted if on average there is an excited (unpaired) spinon sitting between each two 
neighbored holons (as illustrated by Fig. El, where an exchange between a pair of holons can gain a minus sign in 
wavefunction. In other words, the holon condensation must break down when the vortex density n„ is comparable to 
holon density (5, far less than n™^'^ at low doping. Such a consideration provides an estimate of the upper bound for 
as 



T„ determined by Eq. Ij^fil) is plotted in Fig. [Slby filled circles (with E^^a determined by the self-consistent mean- 
field theorji24 of the UPP as ^ 0). Eq. ^ can be also understood based on the "core touching" picture of 
spontaneously excited spinon vortices in the LPP— . Note that the average distance between excited spinons may be 
defined by Is = laj ^Jwn^. Since the characteristic core size of a spinon- vortex is 2ac as shown in Fig. 01 then one 
expects that the "supercurrents" carried by spinon-vortices are totally destroyed when the sample is fully packed by 
the "cores" with 1^ — 2ac which results in Eq. JSHJ. Here it is emphasized that there is no simple mean-field phase 
transition for the LPP due to the peculiar mutual duality gauge structure mediated by A** and A'\ Furthermore, 
27r vortices of ■0/1, independent of spinon-vortices, will be also generally present near as ph — > 0. All of these 
fluctuation effects can affect the crossover from the LPP to UPP in details. 

Similar to the UPP case^i, an external magnetic field will break up some RVB pairs through the Zeeman effect to 
create more excited spinons at a given temperature. By considering the Zeeman effect on the energy spectrum Emcr 
in Eq. H23() . the magnetic field dependence of = Ty{H) can be further obtained from Eq. (|36|l . Or conversely, 
for each T < there is a characteristic field Hy{T) at which the LPP phase is destructed. = H^iT = 0) 
determined this way is shown in Fig. [S] (filled squares). Roughly we find k^T^ 1.5/XB-ff^ according to Fig. [S] or 
Hy/Ty ~ 0.99 (T/K) which is slightly larger than the experimental ratio ^ 0.62 (T/K) determined by the Nernst 
measurement^'' in the LSCO compound. Furthermore, Hy as a function of temperature is further shown in Fig. [7| 
at doping concentration 6 = 0.125 (in the same figure, the critical field Hpq{T) for the UPP is also shown, which 
terminates at T = Tq). 

Superconducting phase. — The superconducting phase transition (the dashed curve in Fig. O is given hy^ 




(36) 



ma 




(37) 



which is determined by the characteristic energy Eg of the low-lying spin excitation — the so-called resonancelike peak 
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FIG. 7: The magnetic field - temperature phase diagram of the pseudogap phases at doping 5 — 0.125. 



energy (c/. Fig. (SJ- Such a transition is non-BCS-like and is represented by a topological transition related to the 
phase coherence of the superconducting order parameter, which will be further discussed in Sec. Ill F. 

In the mixed state below Tc, by including the Zeeman energy according to Ea. iP^ . Eg is reduced to E* — Eg{B) — 
2fiBB such that one can estimate Tc{B) by using a simple relation Tc{B) E*/4,k^. Then in turn one may define 
i7*2 = B{Tc), which is shown in Fig. [3by a dashed curve (here a self-consistent mean-field solution of E* in the LPP, 
in the presence of the Zeeman term, is numerically determined). Note that H*2 so defined will vanish at TJ,, resembling 
the conventional Hc2 in a BCS superconductor. However, since free spinon- vortices are generally present at H > H*2, 
H*2 is a crossover field which no longer has the same meaning as Hc2 in the conventional BCS superconductor. 
Roughly speaking, the Abrikosov magnetic vortices (their total number is proportional to the magnetic field) are 
expected to be present below H*2 where the additional spinon-vortices generated by the Zeeman effect, with both 
signs in the vorticity, are still loosely paired, while the vortices unbinding occurs above H*2 . So H*2 defines a crossover 
between two types of vortex regime. Of course, if the Abrikosov vortex lattice are melt below H*2, the distinction 
between the two vortex liquid regimes should be further blurred at H*2. At T = 0, H*2{0) obtained by setting E* = 
is found to be numerically quite close to i?" at S = 0.125 (Fig. [TJ, although there is no physical reason for them to 
be the same. The numerical result also shows that /iB^^*2(0) — Eg{B — 0,T = 0)/2, or 

h:2{o) (38) 

according to Eq. H37II which is in quantitatively agreement with the critical magnetic field extrapolated from the 
specific heat measurementii in the SC phase by fitting the VolovikA/B term of a d-wave superconductor in the 
presence of a perpendicular magnetic field. 

Therefore, the LPP is distinct from the UPP by an additional charge condensation, and the former is always nested 
inside the latter because, in order to realize the condensation in the charge sector, the excited spinon-vortex density 
should be restricted to avoid too strong the frustrations on the former — each unpaired (excited) spinon always exerts 
destructive phase interference on the charge condensate in a form of carrying a tt flux tube. Here one clearly sees 
how the mutual duality of the gauge structure plays a crucial role in fixing the phase boundary. Finally, a three- 
dimensional phase diagram of the UPP, LPP and SC phase in the parameter of magnetic field, doping concentration, 
and the temperature determined by the above methods is summarized in Fig. |S| 



14 



1- 



1.0 



0.8 



0.6 



0.4 



0.2 




0.0 



FIG. 8: The phase diagram of the pseudogap regimes in the three-dimensional space of magnetic field, doping, and temperature. 



In the absence of the integrity of fermionic quasiparticles (see Sec. Ill F), the transport properties in the LPP are 
mainly associated with the motion of spinon-vortices. This kind of transport is drastically different from the usual 
normal state transport contributed by quasiparticles, which more resembles that in the flux-flow phase of a BCS 
superconductor. But there are some crucial distinctions between the latter and the LPP. In a conventional type II 
superconductor, the vortex core energy is usually very large compared to the temperature T ^ Tc, which means that 
the thermal fluctuations of vortices can be ignored and nearly all the vortices are generated by an external magnetic 
field such that riy = B/(j)Q. On the opposite, in the present LPP, the spinon-vortices as quantum objects can be 
created thermally, which means ^ B/ipo in weak field case. It is those spontaneous vortices that will decide the 
transport phenomenon in such a region. 

Let us start with the charge current given by the London-like equation in Eq. (|30|) . By using the steady current 
condition 



and the definition of the electric field E = —dtAf in the transverse gauge, one can finally obtain the following basic 
relation 



where denotes the density of spinon-vortices and -antivortices with drifting velocity w± along a direction perpen- 
dicular to the electric field in the plane. As illustrated by Fig. Eta), the electric field and the drifting of vortices 
and antivortices must be balanced according to Eq. (|4U|I in order to avoid the system being accelerated indefinitely 
((9f J ^0). So the applied electric field will drive the vortices and antivortices moving along a perpendicular direction 
with opposite velocities: v+ = — v_ if the vortices and antivortices are not polarized by the external magnetic field, 
i.e., = n~ . Equation H4U|) is a basic relation for the spinon- vortex-related transport phenomenon in the LPP. 

The Nernst signal corresponds to the electric field measured along the y-direction when spinon-vortices and - 
antivortices are both driven by a temperature gradient in the same direction along the ^-direction. Such a case is 
shown in Fig. E^b), where the spinon-vortices and -antivortices move along the z-direction with the same velocity 
v+ = v_ = V in a steady state. To have a finite E in terms of Eq. (|40|) . then the vortex density has to be 
polarized by the external magnetic field B =Bz according to the "neutrality" condition H32|l : 



D. Nernst effect 



dtJ = 



(39) 




(40) 




(41) 



such that 



E = zx^o [n, 



— n, 



'V 



) V = B X V 
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FIG. 9: Schematic picture of (a) the flux flowing under an electric field in the LPP, which can lead to a flux flow resistivity as 
well as the spin Hall efi'ect, and (b) the flux flowing under a temperature gradient VT, which must be balanced by an electric 
field E and thus leads to the Nernst effect. 



Define the Nernst signal cat by 



Suppose S0 is the transport entropy carried by a spinon vortex and r/^ is its viscosity such that the drift velocity v" 
can be decided by s^VT = — ?7sV. Then one has 

CN = (43) 



On the other hand, in the absence of the temperature gradient, a charge current can also drive a transverse motion 
of spinon- vortices and -antivortices along opposite directions, i.e., vj- = ±v, such that an electric field is generated 
along the current direction according to Eq. H4U|I . leading to a finite resistivity due to the presence of free vortices, 
which is given by^ 



P = ^^l (44) 

'Is 



where 



n+ + n~ 



denotes the total density of vortices and antivortices. This formula is familiar in the vortex flow regime of a conven- 
tional superconductor. Then, by eliminating rjs which is harder to calculate, one obtains 



= 1^ (45) 

where a^y is the quantity introduced in Ref.^^. It is noted that in the vortex phase of a conventional superconductor, 
one has B = (pon^ = 4>on^ with n~ = such that Eq. H45|l reduces to axy4>o = s,p from which the transport entropy 
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FIG. 10; The quantity aa:y4>o/d, = ejv^o/prf, which is related to the Nernst signal without involving the viscosity coefficient, is 
shown as a function of the magnetic field and temperature (the inset) at a given doping. Here d — 7.7 A is the distance between 
two Cu02 layers. 

of a vortex can be obtained directly (e.g., Ref^). But in the LPP, axy<t'o = s^{B /(pQiiy) < s^, which means the 
transport entropy is underestimated by a^yipo here. 

So far all the derivations are similar to the flux flow region in a conventional superconductor except that the total 
vortex number ny in Eq. (|45|l generally is no longer simply proportional to the applied magnetic field B, which 
remains finite even in the limit B to account for thermally activated spinon-vortices. What really makes the 
Nernst transport unique in the present theory is that the transport entropy here is due to the spin degree of freedom 
associated with a spinon-vortex, instead of a normal core in a conventional BCS superconductor. Let us recall that 
the core size of a spinon-vortex is determined by the "cyclotron" length scale Qc- If one neglects the "configuration" 
entropy associated with the spatial location of a spinon-vortex, then the main "transport entropy" as carried by a 
spinon-vortex is due to its free S = 1/2 moment which gives rise to 

= fcB{ln[2cosh(/3^BS)] -/3^BStanh(/3/iBB)} . (46) 

Here the interaction between spinous in different vortex cores are ignored, which will in general lead to an additional 
suppression of the average entropy per vortex. Under such an approximation, — > 1/2 for T ^ ^gH and ^ 
for T <C hbH. Despite of the saturation of at high temperature, the quantity axy defined in Eq. H45|l decreases 
at high temperature due to the quick proliferation of the total vortex number. The temperature and magnetic field 
dependence of axy4>o/d is shown in Fig. 1101 in which d = 7. 7 A is the distance between two Cu02 layers. The 
magnitude of such a quantity is quite comparable to the experimental datci^^, implying that the transport entropy 
due to the free moment in a spinon-vortex is capable to produce the Nernst signal observed experimentally. However, 
at strong magnetic field or high temperature, axy decays slower than that seen in the real systems. That is because 
here we have only considered single rigid spinon-vortices without taking into account the interaction between spinon- 
vortices as well as the fluctuations of vortices, which become very important as the LPP will eventually collapse at 
Hy and T„ as discussed before. 

Finally, we briefly mention a unique prediction related to the spinon-vortex motion driven by an external electric 
field. As shown in Fig. Ufa), vortices can be driven by an in-plane electric field to move along the transverse direction. 
Since each vortex carries a free moment, if these moments are partially polarized, say, by an external magnetic field 
via the Zeeman effect, then a spin Hall current can be generated along the vortex motion direction. The spin Hall 
conductivity has been recently obtained^ as follows: 
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s f^Xs ( B 

which only depends on the intrinsic properties of the system Uke the uniform spin susceptibihty Xs , with the electron 
g-factor g ~2. It is important to note that the external magnetic field B applied perpendicular to the 2D plane reduces 
the spin rotational symmetry of the system to the conservation of the 5^ component only, satisfying + V • J* = 0. 
Thus the polarized spin current J* is still conserved and remains dissipationless as the current of its carriers - vortices 
is dissipationless in the LPP. In contrast, the charge current remains dissipative with the resistivity given by Eq. H44|l . 
The dissipationless spin Hall effect in the LPP is therefore an important prediction of the phase string model. For 
more details one is referred to Refi^. 




E. Diamagnetism 



In an ordinary superconductor, the diamagnetism can be well described in terms of the GL free energy. In such a 
description, the vortex phase below Tc, as induced by magnetic field, exhibits strong residual diamagnetism, continu- 
ously evolving from a perfect Meissner state and eventually vanishing as the magnetic field H — )■ Hc2 or temperature 
T ^ Tc when the Abrikosov vortices become closely packing ("core touching"). The diamagnetism here is directly 
connected to the energy cost for creating a vortex which carries swirling supercurrents and is subject to long-range 
interactions with other vortices, where the total number of vortices is usually proportional to the applied magnetic 
field by n„ = B/(j)Q. 

The present LPP is a vortex liquid state above Tc, and it would be natural for one to conjecture that a substantial 
diamagnetism still persist. However, like the Nernst effect discussed above, due to the involvement of spontaneous 
vortices as well as the spin degrees of freedom, the theoretical description of the diamagnetism will be also modified 
drastically. In the present of thermal vortices, a weak magnetic field will mainly polarize the numbers of vortices and 
antivortices, instead of creating new vortices, and thus the diamagnetism due to the energy cost of the vortex creation 
will diminish quickly. On the other hand, if there is a dense population of thermal vortices, then the configuration 
entropy will sensitively depend on the imbalance of the numbers of vortices and antivortices, due to the "hard-core" 
like repulsion between two vortices of the same sign of vorticity as discussed before, which will cause an additional 
diamagnetism response. 

Let us note that the total free energy for the LPP is given by 

F^Fh + Fs + I Sv^ (48) 

where B is the magnetic field perpendicular to the plane and Fs is the free energy for the spin part: Fg — 
—f}~^ h\Tr (^e^^^") . The magnetization M — {B — iJ)/47r is determined by 

where denotes the 2D area of the system. 

As discussed before, Eq. (|49|l is composed of two kinds of sources for diamagnetism. The first one is similar to 
that in the conventional vortex phase of a BCS superconductor contributed by supercurrents associated with vortices. 
But the vortex liquid phase above Tc is represented by the presence of a sufficient number of thermally excited free 
vortices that destroy the superconducting phase coherence. These spontaneous free vortices are already present in 
the absence of an external magnetic field perpendicular to the 2D plane, due to the screening of the long-range 
(logarithmic) interaction among vortices [Ea. (|33|l ]. So in the vortex liquid phase, the free energy F, to the leading 
order approximation, should only depend on the total number of vortices and antivortices, each with the same 
renormalized self-energy. On the other hand, when the applied magnetic field is weak enough as compared to the 
total number per unit area of thermally-excited vortices multiplying the flux quantum 0o (i.e., B <^ Uycjio), it can 
only partially polarize the numbers of -t- and — vortices by Eq. (|41|l : An^ — B/(f)Q, while the total vortex number n„ 
remains unchanged. Without an explicit dependence of F on the polarization An^ of vortices and antivortices that 
is proportional to the magnetic field, one would find a vanishing diamagnetic response to the weak magnetic field. 

Another origin of diamagnetism is related to the entropy part of the free energy, as already pointed out above. In 
the following we provide an estimate of this effect as a lower bound of the diamagnetism at weak field limit. Recall 
that we have previously discussed the quantum mechanical nature of the core of a spinon- vortex: the trapped spinon 
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FIG. 11: The diamagnetic susceptibility versus temperature in the LPP. 



does a cyclotron motion within the core according to H^, and its Hilbert space is restricted to the lowest Landau 
level (LLL) at low temperatures (> Tc). The above-mentioned short-range interactions will then dramatically alFect 
the entropy of spinon- vortices: two vortices with the same vorticity will cost a large, unscreened energy if they stay 
in the same real-space (symmetric) cyclotron orbit labelled by m, whereas two with opposite vortices approximately 
do not experience any substantial short-range interaction even when they are in the same real-space orbit. Then we 
can approximately write down the total spinon-vortex configuration entropy as 



ln(-^-l)-.gln(l 



G 



(50) 



where N is the number of lattice sites and = n^A^ is the spinon number with spin index a. G = 6/2x2 = 6 
denotes the degeneracy of the LLL per unit area, including the spin index, and the condition that two same vortices 
cannot occupy the same state labelled by m has been used with assumption that at T > Tc mainly the LLL is involved. 

When the magnetic field is in the range of _B ^ n^4'o, one can reasonably assume SB = {Sn^ — Sn^) = 2(j)oSn^ = 
—2(f)QSn~ since Sn^ — Sn^ + Sn~ ~ 0. Then, the magnetization determined by the entropy Sg is given as follows 



T fdSA T 



n U^Ah, " "(2-^0) {dn+), 

At weak field limit, the rhs of the above expression may be expanded to the first order of B as 

M, = -xSS, or Ms = XdH 

with 

Xci = -T-^ (51) 

where Xd = {2G^n^y^ 4,'^ ■ ^^S- ^2 Xd is shown as a function of the temperature at (5 = 0.125, whose magnitude is 
comparable to the experimental result'^^ at = 1 T. Note that at T = one has = and x^d ~^ such that 
Ms — (I/Att) H, where the full Meissner state is recovered. In reality, one has = below Tc and thus expects to 
see a divergence in Xd there. Near Tf, the correction from the superconducting phase fluctuations may be important 
to result in some nonlinear singular field dependence^^ of Ms which is not included in Eq. (|51|l . Furthermore, Xd 
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reaches a minimum but does not vanish at where ^ S because the rigid spinon-vortices are considered here 
without taking into account of the collapse of spinon-vortices close to T^. Finally, as remarked before, here we have 
only focused on the weak field limit B ^ ny(l)o. With increasing the magnetic field, more spinon-vortices will be 
generated and nonlinear effect at strong magnetic field will emerge which is explored quantitatively elsewhere^. 

F. Superconducting phase transition 

With A'^ 7^ in the LPP, the superconducting state can be finally realized through the phase coherence^! 

^e'*°("")^^0 (52) 

in which the spinon-vortex composites must form vortex-antivortex pairs [cf. Fig. El^b)] in an analog with the KT 
transition in the XY model. 

The corresponding dynamics is decided by Eq. H33|l . It is similar to the XY model, except that A'' introduces tt 
instead 2tt vortices and the vortex cores are attached to spinous which have their own dynamics as governed by Hs 
with an intrinsic quantum length scale ac- Consequently the superconducting phase transition is not of a BCS-type 
second-order one as described by an ordinary GL theory — here iph or more precisely the amplitude y/ph-, is in general 
not expected to vanish at T^, especially in the underdoped regime. It is rather a phase coherence transition involving 
spinon-vortices binding/unbinding which directly affects e'* The detailed renormalization group analysis -^ leads 
to the Tc relation (|37|l which connects the phase coherence temperature with the spin resonancelike energy Eg in a 
quantitative agreement with experimental measurements^^. 

The superconducting phase coherence (|52|l implies that spinous are confined in the bulk where a single spinon- 
vortex excitation costs a logarithmically divergent energy. In this case, the low-lying elementary excitations are S ~ 1 
spin excitations composed of pairs of spinon-vortices. Similar to the LPP, the spin resonancelike mode is basically 
composed of a pair of spinous with Eg = 2Es ^ 6Js with Eg = {Em)uiin denoting the lowest level in the spinon 
spectrum in terms of Hs in Eq. (|21|l . which appears in the dynamic spin susceptibility at Qaf — (t^iT^) as shown in 
Fig. ISJa) in the LPP. However, in the SC phase, the additional correction from Eq. (|33|l can modify such a low-lying 
mode by that for a pair of spinous with different m's, the confinement force from Fh in Eq. H33|l can further add a 
tail in the resonancelike peak from 2Es towards higher energy. 

Another stable elementary excitation in the SC phase is the quasiparticles composed of spinon-holon pairs as dis- 
cussed in Ref.^^ . The d-wave nodal quasiparticle dispersion is schematically shown in Fig. 1121 where the quasiparticle 
spectral function has a sharp lineshape below the spinon characteristic energy scale Eg = Eg/2 as energetically it 
cannot decay into a spinon and a holon in the condensation even locally. Above Eg, however, the composite (spin- 
charge separation) feature can show up in the lineshape of the spectral function, since energetically the quasiparticle 
can decay into a pair of spinon and holon locally without costing much from the logarithmic confining potential. 

Therefore, at low energy and long-distance, there is no more electron fractionalization in contrast to the LPP. A 
single spinon-vortex composite can only emerge as a magnetic vortex core in the mixed state, where it ensures the 
minimal flux quantization at 0o = hc/2e as discussed in Ref.?^. The phase string theory predicts that a free S = 1/2 
moment must appear in the core of a Arbrikosov magnetic vortex, even though a Kondo screening effect due to the 
coupling to the background quasiparticles may complicate the analysis of possible experimental observations of such 
moments. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have systematically explored the physical properties of the lower pseudogap phase in the phase 
string model. In the global phase diagram of the phase string model, the lower pseudogap phase corresponds to the 
regime where the amplitude of the superconducting order parameter becomes finite, but the order parameter itself 
is still short of phase coherence. The presence of free vortices in this phase prevents the system from undergoing 
superconducting until some lower temperature is reached where the vortices and antivortices form bound pairs. So 
the lower pseudogap phase and the superconducting phase are intimately related and the former sets a stage for the 
latter to emerge. 

What makes the present microscopic theory really unique is that the above picture is naturally embedded in a 
spin-charge separation description. That is, for a doped Mott insulator, spin and charge correlations evolve distinctly, 
especially at low doping. Thus, at different energy/temperature scales, the spin and charge degrees of freedom play 
rather different roles. For example, at low doping and high temperature the neutral spins start to form the singlet 
pairing to lower the superexchange energy, which occurs at the expense of the kinetic energy of doped holes and the 
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"composite" lineshape^^*^ 

resonancelike spinon energy 



sharp nodal QP excitation 




FIG. 12: A schematic nodal quasiparticle dispersion in the superconducting phase. The lineshape below the characteristic 
spinon energy -Eg/2 is very sharp as energetically the quasiparticle cannot decay into a spinon and a holon in the condensate. 
Above Eg/2, however, such a decay is energetically allowed locally with the spinon and holon remain loosely confined at large 
distance. 



latter behave quite incoherently in this upper pseudogap phase. So the spin correlations are the dominant driving 
force in this high-temperature pseudogap phase, where the antiferromagnetic fluctuations continuously increase with 
the decrease of the temperature, and would grow into a true long range order at zero temperature if there were no 
intervention from the charge degree of freedom at some lower temperature at finite doping. 

However, the kinetic energy of holes would be highly frustrated in a long-range antiferromagnetically correlated 
spin background. Actually the kinetic energy will eventually win over the long-wavelength spin correlations beyond 
a small finite doping concentration at sufficiently low temperature (Fig. 1). In this regime the trend of growing 
antiferromagnetic fluctuations will get stopped and the system makes a crossover into the lower pseudogap phase 
at low temperature instead. Here the equal-time spin correlation length becomes truncated in the same order of 
magnitude as the average hole- hole distance [c/. Ea. H27l) ]. The corresponding spin gap behavior in the spin channel 
is clearly illustrated by the emergence of a resonancelike peak at energy Eg in the spin dynamic susceptibility at the 
antiferromagnetic wavevector in Fig. Ota) and the suppression of the uniform spin susceptibility at low temperature 
in Fig. Olb). 

The spin liquid behavior in the lower pseudogap phase is thus the consequence of the competition between the 
charge and spin degrees of freedom, in which the kinetic energy of doped holes wins. Though suppressed in density, 
the thermally activated spin excitations across the spin gap can still exert frustration effect on the kinetic energy of 
doped holes as the reminiscence of the competition between the charge and spin sectors, which is the cause for the 
phase disordering. So the superconducting phase can be finally realized only when those spin excitations which bring 
destructive effect on the charge condensate are all excluded out of finite-energy spectra by some form of confinement. 

Technically, in the phase string model represented by the effective Hamiltonians and ©, the spin and charge 
degrees of freedom are mathematically characterized by bosonic spinon and holon fields, respectively, which are coupled 
together by a mutual duality gauge structure. These bare spinous form RVB pairs, and spinon excitations are the 
result of the broken RVB pairs. The lower pseudogap phase corresponds to the condensation of bosonic holons, which 
obviously is connected to the kinetic energy gain in the charge sector. The holon condensation then feedbacks to 
the spin sector to lead to the suppression of low-lying antiferromagnetic fluctuations and the opening of the spin gap 
through the topological gauge structure. The holon condensation is also self-consistently strengthened by this spin 
gap opening effect, because the fractionalized spin excitations-spinons carry fictitious n flux tubes as perceived by 
the holons and thus always introduce destructive phase interference for the latter. 

In the lower pseudogap phase, even though the density of states get suppressed, spinous carrying tt flux tubes still 
can be excited thermally across the energy gap Eg = Eg/2. It is well known that a tt flux tube in a Bose condensate 
will always induce a vortexlike supercurrent response. That is precisely what happens in the lower pseudogap phase 
and each excited spinon is thus locking into a current vortex to form the spinon-vortex composite as shown in Fig. 
^ Therefore the lower pseudogap phase is a vortex liquid/spin liquid phase, which is also known as the spontaneous 
vortex phase since these spinon-vortices can be thermally excited without involving external magnetic field. 

So the central prediction of the present microscopic theory is the existence of spinon-vortices as elementary exci- 
tations in the LPP. A spinon-vortex plays a dual role here. On one hand, it carries a neutral spin-1/2 and indicates 
the spin-charge separation or electron fractionalization in the lower pseudogap phase. On the other hand, it is bound 
to a charge current vortex and contributes to a 27r phase vortex in the superconducting order parameter. The latter 
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property distinguishes the lower pseudogap phase from the upper pseudogap phase, and represents a mid-step to- 
wards the true spin-charge recombination, which is estabhshed in the superconducting phase where spinon-vortices 
are "confined" to form conventional S = 1 spin excitations and to result in coherent quasiparticle excitations by the 
holon-spinon confinement. 

The vortex core of a spinon-vortex is distinctly different from a conventional vortex (like in the BCS supercon- 
ductor). In the latter a competing "normal state" is nucleated at the core, whereas in the former an elementary 
fractionalized object, spinon, naturally serves as a core. These spinon-vortices are actually always present, even in the 
superconducting ground state where all of them are just paired up, instead of "being annihilated". So to create a un- 
paired vortex at low energy/temperature, one only needs to break up an RVB pair without "nucleating" an additional 
"core" . To the leading order of approximation, the energy scale for creating a vortex is the spinon energy gap Eg ~ SJ 
at low doping. In this sense, the spinon-vortices are really "cheap" vortices which can be thermally (spontaneously) 
excited in the lower pseudogap phase. In a similar sense, the lower pseudogap phase is not a "competing" phase of 
essentially different symmetries with regard to the superconducting state. Rather the former is distinguished from the 
latter mainly in the long- wave length limit, much larger than the core size Uc and the RVB pairing size or equal-time 
spin correlation length, where the thermally excited spinon-vortices form vortex-antivortex pairs to realize the 
superconducting phase coherence. 

Consequently the onset of the LPP at Ty or Hy, though a charge behavior, is determined by the spinon excitation 
according to Eq. H36|) . which can be understood based on the "core touching" picture of spinon vortices with the 
"supercurrents" being totally squeezed out. The Nernst effect has been also connected to the spin degrees of freedom as 
the latter provide transport entropy, while the diamagnetism in the low field limit is associated with the configurational 
entropy of spinon-vortices. The overall theoretical results are found to be in the correct magnitude as compared with 
the experimental measurements in the cuprates. Furthermore, a conserved, dissipationless spin Hall effect has been 
also predicated as a unique feature for the LPP. We shall present more detailed and quantitative results based on the 
mutual Chern-Simons gauge theory of the phase string model in a separate publication^. 
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